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Abstract 

We consider ultracold atoms in a two-dimensional optical lattice of the dice geometry in a tight- 
binding regime. The atoms experience a laser-assisted tunneling between the nearest neighbor sites 
of the dice lattice accompanied by the momentum recoil. This allows one to engineer staggered 
synthetic magnetic fluxes over plaquettes, and thus pave a way towards the realization of topologi¬ 
cally nontrivial band structures. In such a lattice the real-valued next-neighbor transitions are not 
needed to reach a topological regime. Yet, such transitions can increase a variety of the obtained 
topological phases. The dice lattice represents a triangular Bravais lattice with a three-site basis 
consisting of a hub site connected to two rim sites. As a consequence, the dice lattice supports 
three energy bands. From this point of view, our model can be interpreted as a generalization of 
the paradigmatic Haldane model which is reproduced if one of the two rim sub-lattices is elim¬ 
inated. We demonstrate that the proposed upgrade of the Haldane model creates a significant 
added value, including an easy access to topological semimetal phases relying only on the nearest 
neighbor coupling, as well as enhanced topological band structures featuring Chern numbers higher 
than one leading to physics beyond the usual quantum Hall effect. The numerical investigation is 
supported and complemented by an analytical scheme based on the study of singularities in the 
Berry connection. 
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I. INTRODUCTION 


Optical lattices have firmly established themselves as a modern and versatile tool to study 
fundamental physics in a clean environment with various physical parameters being under 
experimentalist’s control and often extensively tunable [1-3]. One is typically interested 
in implementing a paradigmatic Hamiltonian that clearly demonstrates a particular phe¬ 
nomenon or an effect. A list of recent successes features, to mention just a few examples, 
realization of the Harper-Hofstadter [4-6] and Haldane models [7], direct observation and 
control of the Dirac points [8], creation of artificial magnetic fluxes via lattice shaking [9] and 
reproduction of models of magnetism [10], engineering of a spin-dependent optical lattice 
resulting from a combination of Raman coupling and radio-frequency magnetic fields [11]. 

In particular, access to topological band structures is of enormous interest [12-14], The 
presence of the topological order is signaled by a non-zero Chern index reflecting a non- 
vanishing integral of the Berry curvature over the entire two-dimensional Brillouin zone. 
A topological band supported by a spatially periodic optical lattice acts as a model of a 
Landau level. The unique band structure consisting of a ladder of Landau levels defines an 
apparent insulator with current-carrying edge states and has traditionally been associated 
with the presence of an external magnetic held. In cold-atom setups, however, the topological 
character becomes an intrinsic property of the band and is not necessarily associated with the 
presence of a physical magnetic held [14, 15]. Synthetic fluxes piercing the lattice plaquettes 
may be imparted by the lattice shaking [9, 14, 16-18], laser-assisted tunneling [14, 19-21] 
or using synthetic dimensions [22], 

Many of the breakthroughs mentioned in the introductory paragraph can be classihed 
as mimicking or reproduction of phenomena known from the condensed matter physics. 
However, significant contributions from cold-atom systems to extending the known physics 
should also be recognized [1, 2, 14, 21, 23]. Perhaps the most obvious examples relate to the 
construction of topological bands with the values of the Chern index greater than one [24- 
30], which is a central topic of the present paper. The properties of such a band is not 
a direct sum of the properties of several Landau levels, and reach beyond the traditional 
(integer or fractional) quantum Hall physics [31, 32], 

Indeed, the study of bands with higher Chern numbers has been particularly relevant in 
connection to the so-called fractional Chern insulators [33-35]. Although many-body inter- 
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actions, which play the central role in these studies, are beyond the scope of the present 
contribution, we stress that many insights into the nature of the fractional topological states 
were obtained from somewhat artificial lattice constructs often involving many layers [25] 
or distant-neighbour hoppings [26, 27, 36]. Ongoing efforts [37-40] are also based on the 
Harper-Hofstadter model that in principle supports subbands of arbitrarily high Chern num¬ 
bers. Here, one also has to defy rather stringent requirements posed by large magnetic unit 
cells, low particle densities, and a large number of subbands implying small topological band 
gaps [40]. In the present paper we focus on exploring the potential offered by relatively sim¬ 
ple and thus more realistic lattice models. We construct a generalization of the Haldane 
model [41-45] by coupling three rather than two triangular sub-lattices. In this way, the 
honeycomb lattice featured in the Haldane model is upgraded to the dice lattice [46-51] 
which supports a three-band model with a clean access to interesting topological configu¬ 
rations, such as bands characterized by the Chern number equal to 2. In the dice-lattice 
model it is just a complex valued nearest-neighbor (NN) coupling that is sufficient to gener¬ 
ate a staggered synthetic magnetic flux and reach nontrivial setups including a topological 
semimetal phase. On the other hand, for spatially periodic hexagonal lattices, non-trivial 
phases can not be reached just by having the complex-valued nearest-neighbor coupling, one 
should add a real-valued next-neighbor coupling [42-44], Note that the dice lattice affected 
by a uniform magnetic flux was used to demonstrate a novel and intriguing mechanism of 
localization of wave packets in Aharonov-Bohm cages [52-54], 

The paper is structured as follows. In Section II, we introduce the lattice geometry and 
derive the 3x3 momentum-space Hamiltonian encapsulating the physics. Then, Section III 
describes the obtained results starting from phases obtained in the presence of NN couplings 
alone and proceeding to more complex configurations requiring next-nearest neighbor (NNN) 
transitions. We conclude with a brief summarizing Section V. 


II. THE MODEL 

A. Lattice geometry 

We consider a dice lattice, which consists of three triangular sub-lattices. One of them 
is called a hub sub-lattice. It is coupled to other two rim sub-lattices, that in turn are not 
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FIG. 1 . (Color online) Left: dice lattice. The blue, green and red sites correspond to three different 
triangular sub-lattices A, B and C. Solid lines show couplings between the sites A and B. Dashed 
lines show couplings between the sites B and C. The primitive lattice vectors are a± and a 2 . Nearest 
sites are connected with the vectors < 5 i, 62 and 63. Right: hexagonal first Brillouin zone of the 
reciprocal lattice defined by the primitive reciprocal lattice vectors b\ and &2- Two inequivalent 
corners are at the points K (red) and K' (blue). 

coupled with each other. Let us denote the hub sub-lattice by B and the rim sub-lattice by 
A and C. The vectors that connect the nearest lattice sites are (Fig. 1): 

<5i — (fix T y/SGy'), S 2 2 V^3<+), £3 (K'x i (l) 

where a is the distance between two such sites. The elementary lattice vectors 

a,i = a(3e x + y/3e y )/2, a 2 = a{ 3e x — y/3e y )/2 (2) 

define a rhombic elementary cell. The set of lattice vectors r n = n\ a \ +77-20-2 (with integers 
? 7 i and ji 2 ) span the hub sub-lattice B (Bravais lattice). The two rim sub-lattices are defined 
in the following way. The first rim sub-lattice A is shifted from the hub sub-lattice B by 
the vector in such way that sub-lattices A and B alone make a honeycomb lattice. The 
second rim sub-lattice C is shifted to the opposite direction by —<5i (see Fig. 1). Let us 
introduce a set of vectors, that span all the lattice sites: 

f’n T sS\. ( 3 ) 
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Here the index s — 0, ±1 labels the three sub-lattices. The sites of the hub sub-lattice 
(s = 0) coincide with the lattice vectors: r n 0 = r n . The sites of the rim sub-lattices A and 
C shifted by ±<5x. i.e. r„ i+1 = r n + Si and r n _i — r n — 8 

It is convenient to introduce an additional lattice vector 0,3 = a 3 — a 2 . The set of the three 
lattice vectors oq (i — 1, 2, 3) together with the opposite vectors —a* connect all next-nearest 
lattice sites, and can be related to 8, as: ai = <5i — 8 3 , a 2 = 8 2 — 8 3 and a 3 = 61 — 8 2 . 

The basic reciprocal lattice vectors 


2,71 /— 2lT /— 

b 1 =—(e x + V3e y ), b 2 = —(e x -V3e 


(4) 


are orthogonal to the lattice vectors, a j • a y = 2n8ij, i,j = 1,2. The first Brillouin zone is 
hexagonal with two inequivalent corners K and K' positioned at K = (2 b 3 + b 2 )/3 and 
K' = (hi + 262)/3. In terms of the Cartesian coordinates these points are given by 


o ht 27 t 

K =—(3e x + V3e y ), K'=—( 3e x ~V3e y ), 


(5) 


as one can see in Fig. 1. 


B. Tight-binding model 


We shall make use of the tight-binding model in which the single particle states |r niS ) 
represent the Wannier-type wave-functions localized at each lattice site r n s , with s — 0, ±1 
being the sub-lattice index. In the language of the second quantization these single-particle 
states read | r njS ) = c^(r„ iS )|vac), where |vac) is the Fock vacuum state, c^(r niS ) and c(r n s ) 
being the creation and annihilation operators of an atom in the corresponding localized 
state. 

The full Hamiltonian of the system consist of three terms, 


// = //,+ II 2 + // 3 


( 6 ) 


The first term H\ describes the laser-assisted tunneling [14, 18-21, 42, 44, 55] of atoms 
between the sites of the hub sub lattice B (s — 0) and its nearest neighboring sites that 
belong to the rim sub-lattices A and C with s — ±1: 

3 


( / f'n)c(v n T s8i) T H. c., 


(7) 


n s=±l 


i=1 
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where J ^ are the coupling amplitudes. Such generalization of dice optical lattice with two 
different hopping parameters and J^ is already discussed in [56]. The laser assisted 
tunneling is accompanied by the transfer of the recoil momentum p s with s = ± 1 , to be 
labelled simply by p± = p±\. In the present situation p + can generally differ from p 
because the transitions between the different sub-lattices can be induced by different lasers. 
Note that the nearest neighbor hopping alone is sufficient to generate flnxes through rhombic 
plaquettes 

$i = ±(p+- P-) ■ a t / 2 , ( 8 ) 

with a* representing a diagonal vector of the plaqnette in question. Yet the magnetic flux 
over the whole hexagonal plaquette remains zero. 

The second term H 2 takes into account the tunneling between the next-nearest neighbor¬ 
ing sites belonging to the same sub-lattice with s = 0 , ± 1 : 

h -2 = J2Y, zY c \ r n,s)c(r ntS + di) + H. c. (9) 

n s=0,±l i= 1 

This term describes the usual (not laser-assisted) hopping transitions between nearest sites 
in each of the three triangular sub-lattices, and J-f' 1 with s — 0 , ±1 are the corresponding 
matrix elements for the tunneling between the atoms belonging to the s-th sub-lattice. 

The third term H 3 describes the energy mismatch for the particles populating the different 
sub-lattices: 

= J2 £ s c j ( r n,s)c{r n}S ). ( 10 ) 

n s=0,±l 

The on-site energies £ s are the diagonal matrix elements of the Hamiltonian in the basis of 
the Wannier states. Without a loss of generality we can take the on-site energy of the hub 
sub-lattice B to be zero: £q = 0. The on-site energies of other rim sub-lattices are to be 
labeled as £±± = £±. 

Since the first term Hi involves complex phase factors that depend on the elementary cell 
number n, the full Hamiltonian H is not translationally invariant. Yet, we will transform 
the annihilation operators according to c(r n 0 ) —> c(r n 0 ) and c(r n s ) —>• c(r n s ) exp(— \p s -r n ) 
with s = ±1, and perform the corresponding transformation for the creation operators. This 
gauge transformation makes the full Hamiltonian ( 6 ) translationally invariant. 

Transition to the reciprocal space is carried out by introducing new operators 

t.HlYpMe*-, c(r n ,,) = ^ w Y / c,(k)e' l ‘ r ~, ( 11 ) 
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together with the Hermitian conjugated creation operators cj(fe). Here IV is a number 
of elementary cells in the quantisation area, and the vectors r n = r n 0 (defined in the 
Subsec. II A) are located at the sites of the hub lattice. In terms of the new operators the 
Hamiltonian (6) splits into its fc-components: 

c + {k) 

H = J2H(k), H(k)= [ct(fe) cj(fc) cl(k)]u(k) c 0 (fc) , (12) 

where 'H(k) is a 3 x 3 matrix: 

e+ + 2J- 2 +) f(k-p+) J (+) g(k -p+/2) 0 

n(k)= j(+)g*(k-p+/ 2) 2 4 0) f(k) J^g(k-p-/ 2) • (13) 

0 J^g*(k-p_/2) s_ + 2jt ) f(k-p.)_ 

Here we also added an extra phase factor to the transformed operators c s (k) —* 
c s (k)e iPs ' sSl / 2 . The functions 

f(k) = E cos (k ■ a t ), g(k ) = e ifc '^ E e- ifc '^ (14) 

i=l i= 1 

entering Eq. (13) are translationally symmetric in the reciprocal space 

f{k + G) = m, g(k + G) = g(k) , (15) 

where G = 71161 + 77-262 is a reciprocal lattice vector, 77-1 and n 2 being integers. Consequently 
the matrix-Hamiltonian 'H(k) is also fully translationally invariant in the reciprocal space 
'H(k) — kL(k + G). Note that Berry curvature in general depends on the choice of Fourier 
transformation (11), while the corresponding Chern number does not [57, 58]. Furthermore, 
the functions f(k) and g{k) obey the following reflection properties 

/(fc) = /(-fc), g{k)=g*(-k). (16) 

All this helps to consider various symmetries of the matrix-Hamiltonian (13). 

III. PHASES OF NON-INTERACTING FERMIONS 

A. Chern numbers and symmetries of the system 

Since the momentum-space Hamiltonian (13) represents a three level system, there are 
three energy bands characterized by energies E n (k ), with n = 1, 2, 3. Each energy band has 
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a Chern number c n to be defined in Eq. (17). We also identify two possible band gaps. The 
first band gap A 12 measures the energy between the first (n = 1) and second (n = 2) bands, 
the second band gap A 23 corresponding to the energy between the second (n = 2) and the 
third (n — 3) bands. 

The Chern number c n for the n-th band is defined in terms of a surface integral of a Berry 
curvature over the first Brillouin zone (FBZ) [14, 59]: 


Cn = d 2 kF n (k). 

2n J fbz 


The Berry curvature F n {k ) can be expressed in terms of the eigenvectors \u n _k) of the 
reciprocal space Hamiltonian (13) as 

(wJ Un ’ k) ) ~ i {w y {Un ' kl ) 

It is well defined as long the eigen-energies E n (k) are not degenerate for any fixed value of 
k. Therefore the Chern number c n can be ascribed to the n-th band if the latter does not 
touch any other bands. If the Fermi energy is situated in a band gap, the Chern number is 
directly related to Hall conductivity due to chiral edge states of the occupied bands [60] via 
a xy = —e^Cn/h [61-63]. For numerical calculation we make use of the discretized version of 
the Berry curvature (18) described in [64], 

For both rim sub-lattices A and C , we set on-site energies of to be symmetrically shifted 
away from the zero point e + = —e_ = e. We also take the tunneling amplitudes to be equal 
= jF) = jA ] — jt°> — ) — j 2 anc [ assume the recoil momenta to be opposite 

p + = —p~ = p for both rim sub-lattices A and C. The choice of opposite recoil momenta 
ensures the maximum flux, because the magnetic flux through a rhombic plaquette given 
by Eq. (8) is proportional to the difference of these vectors. Under these conditions, the 
matrix representation of the fc-space Hamiltonian becomes 

£ + 2 J 2 f{k-p) Jg{k-p/2) 0 

n{k)= Jg*(k — p/2) 2 J 2 f[k) Jg(k+p/2) • (19) 

0 Jg*(k + p/2) -£ + 2 J 2 f{k + p) 

This form of the Hamiltonian exhibits some symmetries. The first symmetry involves in¬ 
version of the on-site energies £—)■—£ followed by the unitary transformation that changes 
the first row with the third one (i.e. interchanges the rim sub-lattices A and C ), as well as 
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9 I \ 

d kJ Un ’ k) 


F n(fc) i 



the momentum inversion k —y — k. Using the reflection properties of the functions / and 
g given by Eq. (16), one arrives at the same Hamiltonian (19). The second symmetry is 
J — y — J, which is a simple gauge transformation. Using these two symmetries we see that 
the change J 2 —» —J 2 gives 'H(k) —» —'H(k). To sum up, all the three mentioned symmetries 
are: (£ — y — £,?f —y 7f), {J —t — J.l-L — y T~C) and ( 1/2 — y ——y —7f). 


B. Numerical analysis 

In this Subsection, we numerically study the Chern phases of non-interacting fermions. In 
order to present dependence of the Chern number on the parameters e, J, J 2 and p we adopt 
a similar presentation of the phase diagram scheme as in [44], We choose the energy unit 
to be the nearest-neighbor tunneling amplitude J. For the recoil momentum p, we express 
the p x component in the units of K x and the component p y in the units of K y , where K is 
one of the FBZ corners, defined in (5). In all the phase diagrams we present the dependence 
of the Chern number c n = c n (p x ,p y ) on the components of the recoil momentum p using 
different colors for each possible values of c n . The areas corresponding to a topologically 
trivial phase with a zero Chern number are shown in green (c n = 0). On the other hand, the 
areas corresponding to non-trivial Chern phases are shown in yellow (c n = 1), red (c n = 2), 
cyan (c n = — 1) and blue (c n = —2). Additionally we display Chern number labels in all the 
presented phase diagrams. 

First we characterize topological properties of the band structure if there is no next- 
nearest neighbor coupling (J 2 = 0). In the Fig. 2 we show the Chern number phase diagrams 
for £ = ./. One can identify regions where the Chern numbers are {ci,c 2 ,cs} = {0,0,0}, 
{—1,2, — 1} and {1, —2,1}. In the first type of the regions (green color) we have topologically 
trivial regions. In other regions there are non-zero Chern numbers with band gaps Ai 2 = 
A 2 3 = 0. Analysis of the band structure in these regions shows that the bands do not overlap 
and touch indirectly. Thus by filling the first one or the first two bands we arrive at semi- 
metallic phase with non-zero Hall conductivity. The typical spectrum of such non-trivial 
semi-metallic case is presented in Fig. 6. The size of the non-trivial regions in the p-plane 
depends on the mismatch £ of the on-site energies of A and C sublattices. By increasing £ 
from zero these regions immediately appear around the points p = K and become larger in 
size. For about £ ~ J these regions have the largest area as presented in the Fig. 2 for e = J. 
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FIG. 2. (Color online) Chern number dependence on the recoil momentum p in the case e = J 
and J 2 = 0. In the left panel we present the phase diagram of the lowest band Chern number c\. 
In the right panel we show corresponding phase diagram for the middle band. Since the sum of 
Chern numbers over all three bands is zero, the third band gives the same phase diagram as the 
first one (ci = C3). The green regions correspond to the Chern number zero. The yellow, red, cyan 
and blue regions correspond to the Chern numbers 1, 2, -1 and -2 respectively. Nonzero Chern 
numbers are also displayed as labels. The hexagon represents the FBZ in the p-plane. 

For even larger values of e the non-trivial regions shrink back to the points K and finally we 
are left only with the trivial phase {0, 0, 0} everywhere. The analytical treatment, presented 
in the section IV gives the value of £ = for which the semi-metal regions completely 
disappear. For J 2 = 0 there are no other types of phases than the trivial and semi-metallic 
discussed above. Nonzero band gaps appear only in the regions of trivial phase. 

For the case J 2 = 0, the change p —)• p + G, where G is the reciprocal lattice vector, 
corresponds to a gauge transformation. Thus there is a symmetry (p — > p + G, T-L —> T-L). In 
the phase diagram (Fig. 2) we also show the FBZ in the p-plane, which is a hexagon with 
two inequivalent corners positioned at the points K and K'. 

Now let us analyze effects of the non-zero next-nearest neighbor coupling. For this we 
set J 2 = 0.3J and £ = 2 J. The phase diagrams of the Chern numbers are presented in the 
Fig. 3. We can see regions with the Chern numbers corresponding to trivial phases {0, 0, 0} 
and phases {0, ±1, =pl} and {±1, 0, =Fl}- In the latter two types of regions we can find points 
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FIG. 3. (Color online) Chern number dependence on the recoil momentum p in the case e = 2 J 
and J 2 = 0.3J. In the left panel we present the Chern number ci of the lowest band, while in 
the right panel we show the Chern number C 2 of the middle band. For the third band (not shown 
here) we have C 3 = — (ci + C 2 ). The green, yellow, red, cyan and blue regions correspond to the 
Chern numbers 0, 1, 2, -1 and -2 respectively. Nonzero Chern numbers are also displayed as labels. 
A smaller hexagon shows the FBZ corresponding to the case J 2 = 0. Since the introduction of 
non-zero J 2 changes the periodicity of the p-dependence, we also show a bigger hexagon, which is 
now the FBZ in the p-plane. 

corresponding to non-zero band gaps A 12 > 0 and/or A23 > 0 (Fig. 4). This shows that 
there exist topological Chern insulating phases. For example at the point p = K, we have 
the Chern numbers {0, —1,1}, the band gap between the middle and highest bands being 
A 23 ~ 0.26J. Band widths in this case are about 3 J. By positioning the Fermi energy 
in the gap between the second and third bands one arrives at the Chern insulating phase. 
Another interesting point is p = 2 K, which gives the Chern numbers { — 1,0,1}, the band 
gaps A 12 ~ 1.55J and A 23 ~ 0.54J and the band widths of about 2 J. The bottom and 
top bands have non-zero Chern numbers, while it is zero for the middle band. Depending 
on the filling there are two types of topologically non-trivial phases. If the Fermi energy 
is positioned in one of the band gaps, we get a topological insulating phase. If the Fermi 
energy is situated within a band, the band is partially filled and supports the Chern metal 
phase. The discussed types of Chern number distributions over the bands are typical when 
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FIG. 4. Dependence of the band gap on the recoil momentum p in the case where e = 2 J and 
J 2 = 0.3J. In the left panel we present the band gap A 12 between the first and second bands. In 
the right panel we show the band gap A 23 between the second and third bands. 

J -2 is non-zero and smaller than J and e. 

In the case of non-zero NNN coupling J 2 the translation symmetry in the recoil momentum 
p is smaller than in the case of zero NNN couplings: one has to shift the momentum by 2 G 
rather than G. In the phase diagram presented in the Fig. 3 we show this by extending the 
FBZ, which is now a bigger hexagon. 

There are more types of Chern phases when the coupling J 2 is larger than in the previous 
discussion and comparable to the on-site energy e. For £ = J 2 = 0.5 J we find insulat¬ 
ing phases with Chern numbers (±1, ± 1 , 4 = 2 } and metallic phases with Chern numbers 
{±2,0, 4 = 2 } (Fig. 5). For example in the point p = 2 K we get Chern numbers C\ = C 2 = — 1 
and C 3 = 2 with band gaps A 12 ~ 0.61 J and A 23 ~ 0.54J. The width of the lower two bands 
are around 3 J, while the band width of the highest band is about 1.5J. Another interesting 
point is p = (2 K x , K y ) where the Chern numbers are cy = —2, C 2 = 0 and C 3 = 2 (white 
point in the Fig. 5). The bulk spectrum in this point is given in the Fig. 6 . Note that there 
is a gap A 13 « 1.35J between the lowest and highest bands. In this gap there is a middle 
band with a zero Chern number. By setting the Fermi energy in this gap one gets the Chern 
metallic phase with the Chern number ci = —2. 

To summarize the numeric analysis for J 2 7 ^ 0, the typical nontrivial Chern num- 
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FIG. 5. (Color online) Chern number dependence on the recoil momentum in the case e = 0.5J 
and J2 = 0.5 J. Left : Chern number ci of the lowest band. Right : Chern number C2 of the middle 
band. The color scheme and labeling are the same as in the figs. 2 and 3. The white point is 
p = (2 K x ,K y ) where the Chern numbers are ci = —2, C 2 = 0 and C 3 (see the spectrum in the 
Fig. 6 ). 

ber distributions over the bands are {0, ±1, q=l}, {±1,=f1,0}, {=hl, 0, =Fl}, {±1,±1,=f2}, 
{±2, =Fl, =Fl} and {±2,0 ,=f 2}. One can also find the case {±1, =p2, ±1}, which is typical 
for J 2 = 0. For smaller J 2 compared to J and e, one usually gets Chern numbers up to 
1 in modulus. For larger Chern numbers (up to 2 in modulus), one needs to make the 
NNN-hopping J 2 be comparable to the on-site energy mismatch e. 


IV. ANALYTICAL CHERN NUMBER CALCULATION 

Analytic Chern number calculation is based on integration of a Berry connection around 
each singularity point. The Berry connection of the n-th band is defined as [14, 59] 

-A n (fc) i(Uk,n | Vfc , (20) 

where \uk, n ) denotes the n-th eigenvector of the matrix (13). One can express the Berry 
curvature (18) as the z component of the curl B n = V x A n , namely F n {k ) = e z ■ B n . Using 
the Stoke’s theorem we change the integral featured in Eq.(17) over the FBZ to a contour 
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FIG. 6 . (Color online) Bulk lattice spectrum projected along k y for a number of different k x values 
in the range — K x <k x < K x . Left-, the spectrum for the recoil momentum p = K in the absence 
of the NNN-coupling (J 2 = 0) and for e = J corresponding to the parameters used in the Fig. 2. 
In that case there is no energy gap in the spectrum, but different energy bands do not directly 
touch each other. A topological semi-metal phase is formed if the atoms £11 the first energy band 
or the first two bands. Right: the spectrum for the recoil momentum p = (2 K x ,I\ y ) in the case 
where e = 0.5J and J 2 = 0.5J corresponding to the phase diagram shown in Fig. 5. Now there are 
two bands with non-zero Chern numbers ±2 separated by a quasi-gap A 13 ~ 1.35J containing a 
middle band with a zero Chern number. 


integral around the FBZ, 
1 

27 T 


f d 2 kF n (k) —tv—/ 

JFBZ 27T J FI 


d k ■ A, 


£ 


dfc ■ A r . 


singul 


FBZ 27T 

where the last term excludes any contribution due to unphysical gauge-dependent singular 
points of the Berry connection [44, 65, 66]. Since the fc-space Hamiltonian H(k), given by 
Eq. (13) or (19), and its eigenstates are periodic in the FBZ, A n is also periodic. Thus 
the contour integral around the FBZ (the first term on the r.li.s. of the above equation) 
is zero. Consequently the Chern number (17) can be calculated by integrating A n around 
each excluded singular point [44]: 


= —£ 
2tt z ^ 


dk ■ A r 


singul 


( 21 ) 


where the sum is over all singular points in the FBZ. 
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Let us summarize our analytical results, details being presented in the Appendix. For the 
case where the recoil momentum coincides with the inverse lattice vector (p = G) we always 
have trivial phase with all three Chern numbers equal to zero. For the semi-metal case 
(Fig.2) with no NNN-hopping and p = K we find two phases, depending on the mismatch £ 
of the on-site energies. If £ < Eq = we get Chern semi-metal phase with Chern numbers 

{1,-2,1}. Ife > £q, we get a trivial phase {0, 0, 0}. In this way at larger mismatch between 
the on-site energies the topological phenomena disappear. This is in agreement with the 
numerical calculation presented in the previous Section. 

It is possible to apply this method for other values of the recoil momenta p and for 
a general non-symmetric case with the NNN-hoppings. In such calculations one needs to 
diagonalize the matrices of the size at most 2x2. Yet generally ordering of the eigenvalues 
might be a quite involved task, especially if they depend on more than one parameter. 


V. CONCLUDING REMARKS 

In conclusion, we have considered a two-dimensional dice lattice operating in a tight- 
binding regime. The laser-assisted nearest neighbor transitions are accompanied by the 
momentum recoils. This allows one to engineer staggered synthetic magnetic fluxes and 
thus facilitates realization of topologically nontrivial band structures. Real valued next 
nearest neighbor transitions - although not necessary in principle to reach the topological 
regime - may also be present and contribute to the richness of the obtained topological 
phases. The considered dice lattice represents a triangular Bravais lattice with a three-site 
basis consisting of a hub site connected to two rim sites, providing three energy bands. Thus 
our model can be interpreted as a generalization of the paradigmatic Haldane model which 
is reproduced if one of the two rim sub-lattices is eliminated. We have demonstrated that 
the proposed upgrade of the Haldane model creates a significant added value such as (i) an 
easy access to topological semimetal phases relying on only the nearest neighbor coupling 
and (ii) enhanced topological band structures featuring Chern numbers higher than one 
and thus providing access to physics beyond the usual quantum Hall effect. The numerical 
analysis have been supported by an analytical scheme based on the study of singularities in 
the Berry connection. 
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Appendix A: Details on analytical Chern number calculation 


1. Momentum space Hamiltonian and its eigenstates 

Let us establish a general structure of the eigenstates for the matrix Hamiltonian 'H(k), 
Eq.(19). For this we introduce a basis of our three-level system |s), with s = 0, ±1, and 
rewrite the matrix Hamiltonian in the state-vector notation as 

n(k)= £ \s)d s (k)(s\+ £ (| S )^(fe)e i -W(0|+H.c.), (Al) 

s=0,±l s=±l 

where d s (k) stands for the diagonal matrix elements: 

d s (k) = se + 2J 2 f(k — sp). (A2) 

The off-diagonal matrix elements 

Jg(k T p/‘2)=g±(k)e ia ^ k) (A3) 

have been represented in terms of their amplitudes g±\{k) = g±(k ) and phases a±\(k) = 
a±(k). 

Since there is no coupling between the A and C sub-lattices, one can perform a k- 
dependent unitary transformation eliminating the phase factors 

|s) —> |s, fe) = |s)e lsas( - fc \ s = ±1, 

and leave the basis vector |0) unchanged (|0) = |0,fc)). In the new basis the Hamiltonian 
(Al) is characterized by real and symmetric matrix elements. Its eigenvectors can be cast 
in terms of these vectors with real coefficients C n:S (k): 

K,n>= £ C n>s (k)\s,k) = £ | S >C£(fc)e is ^ fc ), (A4) 

s=0,=tl s=0,±l 

Combining Eqs. (20) and (A4), one arrives at the following expression for the Berry connec¬ 
tion 

A n (k) = - £ sCl s (k)V k a s (k). (A5) 

s=±l 

This result together with Eq. (21) will be subsequently used in Ending the Chern numbers. 
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2. Determination of the Chern numbers: General 


To determine the Chern number given by (21), one needs to study a behavior of the 
vector potential at its singular points. Singularities of the vector potential can emerge at 
the points where the phase of the coupling matrix element g±(k)e ia± ^ given by Eq. (A3) is 
undefined. This happens if the function g(k —p±/ 2) goes to zero. The function g{k) given 
by Eq. (14) is zero at the corners of the FBZ, namely at two inequivalent points K and K'. 
Thus there are two pairs of points 

K ± = ±p/2 + K , K' ± — ±p/2 + K' . (A6) 

at which the function g(k =F p/2) goes to zero and its phase a±(k) is undefined. Let us 
determine the coupling matrix element g±(k)e ia± ^ in a vicinity of these points. Combining 
Eqs. (14) and (A3), the amplitude and phase of the coupling element reads up to the first- 
order in the displacement vector q , i.e. for qa <C 1 with q — \q\: 

g±(K± + q) « \qaJ , a±{K± + q) « § - <p, (A7) 

g±( K ± + q) ~ IqaJ, a±{K'± + q) « -f + if, (A8) 

where ip is a phase of the complex number q x + i q y = qe lv> . Integrating over a small circle 
centered at q = 0 surrounding each singular point of the phase, one finds: 



d q ■ V q a ± (K± + q) 


—2vr, 



d q ■ \7 q a±{K' ± + q) 


2tt, 


where the signs are different due to the opposite phases in Eqs. (A7) and (A8). These 
equations together with Eqs. (21) and (A5) provide the following result for the Chern number 


Cn = Y, s \Cl,,(K,) - Cl,(K[) 1 , (A9) 

S=±l 

with K± i = K± and K' ±1 = K’ ± . Therefore to find the Chern number one needs to 
determine the coefficients C n>s entering the state-vector at the points K± and K' ± . If 
C%± = 1, the corresponding singular point contributes to the Chern number of the n-tli 
band. In the following we shall analyze two different situations. 
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3. Determination of the Chern numbers: Specific cases 


Since the Hamiltonian TL(k) given by Eq. (19) or (Al) has a symmetry (e —>■ —e, Ti —» TL), 
we consider only the case where £ > 0. 

a. The case where p = G 

Suppose first that the difference in the recoil momenta coincides with the inverse lattice 
vector p = G. In that case the coupling completely vanishes both for k = K± and also for 
k = K±. At these points g{k — p/2) = g(k + p/2) = 0, so all the states s) (s = 0, ±1) 
are decoupled, and thus the eigenstates are the bare states s). The corresponding eigen- 
energies of the matrix Hamiltonian 'H(k), Eq.(Al), coincide with its diagonal elements 
d s (k) for k = K± and k = K' ± . Since p = G, one has f(k — p) — f(k) = f(k + p), 
giving d s {k ) = se + 2J-2f(k). Therefore the eigenstates are ordered in the same manner 
d+i{k) > d 0 (k) > d_i(fc) both for k = K± and also k = K±, giving C'^ S (K S ) = C/ 1S (K' S ) 
with s = ±1. As a result, the Chern number given by Eq. (A9) is identically equal to zero, 
and the system does not exhibit any topologically non-trivial phases. This is because for 
p = G the flux over the rhombic plaquettes <fh = ±p ■ di is zero (modulo 2n), and there is 
no breaking of the time-reversal symmetry. 


b. The case where p = K 

As another illustration we pick the recoil momentum p = K and take J -2 = 0. In that 
case the Chern numbers have been numerically found to be c± = 1, C 2 = —2 and C 3 = 1, 
see Fig. 2. By taking p = K the phase singularities of the coupling elements g(k T p/2) 
emerge at the points K± = ±K/2 + K and K' ± = ±K/2 + K ', as one can see in Fig. 7. 
Furthermore, the point k = K + is equivalent to the point k = K' . For the latter two 
points we have g(k — p/2) = g{k + p/2) = 0, so there are no coupling matrix elements. 
Since J 2 = 0, the Hamiltonian (Al) at these points is simply 

H{K+) = H{K'_) = £ s\s)(s\, (A10) 

s=±l 

so the diagonal energies entering the Hamiltonian (Al) are d s {k) = se. 
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FIG. 7. (Color online) The phase singularity points K± and K± of the coupling matrix elements 
g(k =F p/2) for p = K. The points K + and K'_ are equivalent. They are shown by red dots 
connected with a double arrow. 

Eigenvalues, ordered from the lowest to the highest, are E l (K + ) = E^K'_) = -e, 
E 2 (K + ) = E 2 (K '_) = 0 and E 3 (K + ) = E 3 (K'_) = e. There is no degeneracy for e > 0 and 
the coefficients C n ^ + (K + ) and C n -(K '_) do not change if one increases s. The only non-zero 
coefficients contributing to the Chern numbers read 

C 3t+ (K+) = C lt .(K'_) = 1. (All) 

For the point k = K _ the non-diagonal matrix elements of (Al) are Jg(k +p/ 2) = 0 
and Jg(k — p/2) = 3 J. Similarly for the point k = K' + these elements are Jg(k — p/2) = 0 
and Jg(k + p/2) = 3 J. Thus the Hamiltonian (Al) at these points is 

H(K-) = £ E„±i s|s)(s| + 3J(|0){+| + |+>{0|), (A 12 ) 

H(K) =£E,^i£|s)(s|+3J(|0)(-|+|-)<0|). (A13) 

Eigenvalues of the H{KJ) are E^ Q \KJ) = — e and E^\KJ) — \ (e ± yj e 2 + 36J 2 ), and 
whose of H(K' + ) are E^\K' + ) = £ and E^(K' + ) = \ e ± yje 2 + 36J 2 ). They are 
plotted in Fig. 8. For e = = £o there are degeneracies E^°\K_) = E^~\K_) = — £ 0 

and E^\K' + ) = E^ + \K' + )=£:o- The eigenvalues change their order at the crossing point 
e = £(u as one can see in Fig. 8. 

Let us first consider the case 0 < £ < £q. The eigenvalues of H(KJ) are in the in¬ 
creasing order: Ei(K_) = E^(K_), E 2 {K _) = E^(K_) and E 3 (K _) = E^+\K_). On 
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FIG. 8. Dependence of eigenvalues of the Hamiltonian H(k) on the on-site energy e for p = K in 
the absence of the next-nearest neighbour coupling. The eigenvalue E ^ is plotted in red dashed 
lines to distinguish it from the other eigenvalues IfU) Left : eigenvalues at the point k = K 
Right: eigenvalues at the point k = K' + . The eigenvalue crossing point e = J = £q corresponds 
to a transition from a topological semimetal phase on the left to a trivial phase on the right. 

the other hand, coefficients required for the Chern number calculation are C\-(KJ) = 0, 
C 2 -(K_) = 1 and C 3 -{KJ) = 0. Similarly H(K ' + ) gives the eigenvalues E\{K' + ) = 
E(~\K' + ), E 2 (K' + ) = EW(K' + ) and E 3 (K’ + ) = E( + \K' + ) and the coefficients C h+ (K' + ) = 
0, C- 2 : +(K' + ) = 1 and C 3t+ (K' + ) = 0. Combining this result together with (All) we col¬ 
lect four non-zero coefficients: C 3t+ (K + ), C\-{K'_), C 2 -(K_) and C' 2 ,+ (K' + ). Substituting 
them into Eq. (A9), we get the Chern numbers for each energy band 

d =Cl_(KL) = 1, 

c 2 = -Cl + (K' + ) - Cl_(K_) = -2, 
c 3 = Cl + (K + ) = 1. 

This result agrees with the numerical analysis presented in Fig. 2. 

Now let us consider the case e > Eq. From the Fig. 8 we see that the eigenvalues 
are reordered as E\ (K_) —* E 2 (K_), E 2 (K' + ) —> E 3 (K' + ), so non-zero coefficients are 
C 3j+ (K + ), and C 3}+ (K' + ). Using Eq. (A9), one can see that the Chern 


(A14) 

(A15) 

(A16) 
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numbers of all bands are now zero: 


ci = —Cf_(-RT_) + Cl_(K'_) = 0, (A17) 

c 2 = 0, (A18) 

c 3 =Cl + {K + )-Cl + {K’ + ) = 0. (A19) 

Thus there is a topological phase transition at £ = J corresponding to the eigenvalue 
crossing in Fig. 8. 
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